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We consider multi-qubit systems and relate quantitatively the problems of generating cluster states 
with high value of concurrence of assistance, and that of generating states with maximal bipartite 
entanglement. We prové an upper bound for the concurrence of assistance. We consider dynamics of 
spin-1/2 systems that model qubits, with different couplings and possible presence of magnètic field 
to investigate the appearance of the discussed entanglement properties. We find that states with 
maximal bipartite entanglement can be generated by an XY Hamiltonian, and their generation can 
be controlled by the initial state of one of the spins. The same Hamiltonian is capable of creating 
states with high concurrence of assistance with suitably chosen initial state. We show that the 
production of graph states using the Ising Hamiltonian is controllable via a single-qubit rotation 
of one spin-1/2 subsystem in the initial multi-qubit state. We shown that the property of Ising 
dynamics to convert a product state basis into a special maximally entangled basis is temporally 
enhanced by the application of a suitable magnètic field. Similar basis transformations are found to 
be feasible in the case of isotropic XY couplings with magnètic field. 

PACS numbers: 03.67.-a, 03.67.Mn, 64.60.Cn 



I. INTRODUCTION 

One of the key features of a physical system for quan- 
tum information processing (QIP) is quantum entangle- 
ment. The problem of entanglement of multipartite sys- 
tems is far from being completely understood, and it has 
numerous interesting aspects. 

One of the possible approaches to multipartite entan- 
glement is to search for quantum states with prescribed 
bipartite entanglement properties Q-Hl-Q This is a non- 
trivial task as there exist limitations on the bipartite en- 
tanglement in multipartite systems, which wcrc quanti- 
fied by Coffmann, Kundu and Wootters 3 . In a pioneer- 
ing work, O'Connor and Wootters Q have considered a 
system of quantum bits, and have searched for an entan- 
gled state of these with maximal bipartite entanglement. 
This state appears to be the ground state of the antiferro- 
magnetic Ising model, the spins representing the qubits. 
This illustrates the relation between states of maximal bi- 
partite entanglement and the spin couplings known from 
statistical physics. We will refer to this approach as the 
question of direct bipartite entanglement, as the relevant 
quantity is the bipartite entanglement present in the sys- 
tem as it is. 

Another approach to the problem of multipartite en- 



*On leave from: Research Institute for Sòlid State Physics and 
Òptics of the Hungarian Academy of Sciences, Budapest 



tanglement is related to cluster and graph Q states. 
These are genuine multipartite entangled states, which 
can be projected onto a maximally entangled state of 
any chosen two spins by a von Neumann measurement 
on the others. Such states arise dynamically in a system 
of spins with pairwise Ising couplings. They constitute 
the fundamental entangled resource for one-way quantum 
computers M, • It is an interesting property of the Ising 
dynamics in this case, that it transforms a whole basis 
of product states into a basis which consists of cluster or 
graph states. In this way a basis transformation from a 
product state basis to a special - in a sense maximally - 
entangled basis is realized. 

These states are the starting points for the second ap- 
proach, the bipartite entanglement in multipartite sys- 
tems available via assistive measurements on all but two 
subsystems. The two key concepts in its quantitative 
description are entan glem ent of assistance [I(J (or con- 
currence of assistance [Hj , quantifying the entanglement 
available via assistive measurements, and localizable en- 
tanglement [Ï3 . IÏ3j . The computational feasibility of 
concurrence of assistance for a pair of qubits makes the 
quantitative study of a part of this question feasible. 

One of our aims is to relate the above two approaches. 
We will show that the optimizations of direct and mea- 
surement assisted bipartite entanglement are indeed re- 
lated. Our other task is to study these genèric features 
in actual spin systems, as such systems do appear quite 
naturally in this context. 

Coupled spin systems have attracted a vast amount 
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of research interest in the quantum information commu- 
nity recently. The couplings studied in statistical physics 
allow for performing certain tasks in QIP such as e.g. 
quantum state transfer 0, 0, ^| , realization of quan- 
tum gates 0, ^|, and quantum cloning 0. As the 
systems of coupled spins are appropriate models for sòlid 
state systems, and also for quantum states in optical lat- 
tices in certain cases [2(|, they bear actual practical rel- 
evance. 

In the second part of this paper we focus on dynam- 
ical generation of entanglement. We consider a system 
initially in a pure product state, and investigate the en- 
tanglement of the states of the system throughout the 
cvolution. The "prototype" of such entanglement gen- 
eration is that of cluster and graph states. The vari- 
ous aspects of the dynamical behavior of entanglement 
in spin systems has been considered by several authors 
recently [U 

In addition to intcrpolating between the two ap- 
proaches to bipartite entanglement in multipartite sys- 
tems, we consider the possibility of controlling the pro- 
cess through the initial state of the system. We address 
the following question. Is it possible to dynamically gen- 
erate states with optimal direct bipartite entanglement? 
We find a positive answer, and also that the same cou- 
plings are capable of producing states with high bipar- 
tite entanglement available via measurements, if a diffcr- 
ent initial state is chosen. Our main tool of describing 
measurement assisted bipartite entanglement will be con- 
currence of assistance. We will examine the possibility 
of controlling the behavior of this entanglement genera- 
tion by the initial state of the system. This is analogous 
to the control of quantum operations in programmable 
quantum circuits 0, 0, 0, 0] ■ Finally we show that 
a suitably chosen magnètic field can enable couplings 
different from Ising to create whole entangled bases re- 
sembling those of cluster states regarding concurrence of 
assistance. (Note that the generation of cluster states 
with non- Ising couplings was considered very recently in 
Ref. 0]) In addition, the application of magnètic field 
in the case of Ising couplings can temporally enhance the 
presence of high pairwise concurrence of assistance. 

As we are mainly interested in illustrating genèric fea- 
tures and certain examples of entanglement behavior, a 
part of our results concerning actual spin systems is sim- 
ply computed by numerical diagonalization of the appro- 
priate Hamiltonians, even though we present some an- 
alytical considerations where we find them appropriate. 
Thus some of our considerations are limited to an order 
of 10 spins, even though according to the numerical expe- 
rience, they seem to be scalable. This number coincides 
with that of the quantum bits expected to be available in 
quantum computers in the near future. As the realization 
of the discussed couplings is not necessarily restricted to 
spins, our results may become directly applicable in such 
systems. We consider two topologies of the pairwise in- 
teractions: a ring where each spin interacts with its two 
neighbors, and also the star topology where the inter- 



action is mediated by a central spin interacting with all 
the others. This was found interesting from the point 
of view of entanglement distribution [32] and also from 
other aspects of its dynamics 0] recently. 

The paper is organized as follows: in the introductory 
Section m we briefly review the entanglement measures 
we use in the following. Section II I II is devoted to the 
review of the dynamical generation of cluster and graph 
states in spin systems, which is the background of the 
second part of the paper. In Section ITVI we present two 
interesting properties of concurrence of assistance, which 
relates the two above mentioned approaches to bipartite 
entanglement in multipartite systems, and will be useful 
in the following. In Section the controlled generation 
of specific entangled states is addressed. Section IVII is 
devoted to the enhanced generation of certain entangled 
bases with the help of magnètic field. Section IVIII sum- 
marizes our results. 



II. ENTANGLEMENT MEASURES 

In this Section we give an overview in a nutshell of the 
entanglement measures and related quantities that will 
be used throughout this paper. 

a. One-tangle. For a bipartite system AA (A being 
a qubit, A being the rest of the system) in the pure state 
|Í") A ^, the one-tangle of either of the subsystems 

T(|*)^)=4det( eA ) (1) 

(where qa = tr A l 1 ! 7 )^ (^1), is a measure of entangle- 
ment. It quantifies the entanglement between the qubit 
A and the rest of the system, including all multipartite 
entanglement between qubit A and the sets all the sub- 
systems in A. 

Although there is an extension of one-tangle to mixed 
states, it is not computationally feasible except for the 
case of 2 qubits, in which case one-tangle is equal to 
the square of concurrence. This justifies the following 
interpretation: the square root of one-tangle is the con- 
currence of such a two-qubit system in a pure state, for 
which the density matrix of one of the qubits is equal to 
that of qubit A. This means, it would be the concurrence 
itself if the subsystem A were also a qubit. 

b. Concurrence. Having a bipartite system in a 
mixed state, a way of dcfming their entanglement is to 
consider the average entanglement of all the pure state 
decompositions of the state. This quantity is termed as 
the entanglement of formation: 

E{q) = win^piEd^i)), so that = g. 

i 1 

This is a kind of generalization of the entanglement de- 
fined in Eq. (JIJ ■ Its additivity is one of the most inter- 
esting open qüestions of QIT. 

The definition of entanglement of formation supports 
the following interpretation: imagine that the bipartite 
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system as a whole is a subsystem of a large system. En- 
tanglcmcnt of formation measures the bipartite entangle- 
ment available on average if everything but the bipartite 
subsystem is simply dropped. 

If the system in argument consists of two qubits, thcre 
is a closed form for entanglement of formation found by 
Wootters [35|. This consideration includes another en- 
tanglement measure. 

Given the two-qubit density matrix g, one calculates 
the matrix 



q= (a (y) ®a (y) )g*{a (y) ®a (y] ), 



(3) 



where * stands for complex conjugation in the product- 
state basis, p describes a very unphysical state for an 
cntangled state, while it is a density matrix for product 
states. 

In the next step one calculates the eigenvalues Ai (i — 
1 ... 4) of the Hermitian matrix 



R = \ \/Q(hfQi 



(4) 



which are in fact square roots of the eigenvalues of the 
non-Hermitian matrix 

i?2 = gg. (5) 

Concurrence is then defined as 

C(g) = max(0, Ai - A 2 - A 3 - A 4 ), (6) 

where the eigenvalues are put into a decreasing order. 
Entanglement of formation is a monotonously increasing 
function of concurrence: 



E(g) = h 



1 + VT= C{gf 



h(x) :— —x\og 2 (x) — (1 — x) log 2 (l — x). (7) 

Thus concurrence can be used as an entanglement mea- 
sure on its own right. 

In multipartite systems the one-tangle and concurrence 
are linked by the Coffmann-Kundu- Wootters inequalities 



(8) 



Ijtk 



which have be proven initially for three qubits in a pure 
state and certain classes of multi-qubit states. For a long 
time they were conjectured to be true in general. This 
conjecture was very recently proven |3jj. These inequal- 
ities set limitations to the bipartite entanglement that 
can be present in a multipartite system. 

c. Concurrence of assistance. Consider again a bi- 
partite system described by the density operator g. One 
can follow a route complementary to that in case of 
entanglement of formation and ask what is the maxi- 
mum average entanglement available amongst the pure 



state realizations, termed as the entanglement of assis- 
tance 35 1: 



i 

so that I*,) = g, 



(9) 



c.f. Eq. @. 

Interpreting again the bipartite system as a subsystem 
of a larger system, one can consider that the whole sys- 
tem is in a pure state, that is, we have a purification of 
g at hand. In this case entanglement of assistance de- 
scribes the maximum entanglement available on average 
in the bipartite system, when a collaborating third party, 
instead of omitting the rest of the system as in the case 
of entanglement of formation, makes optimal von Neu- 
mann measurements on it. Although entanglement of 
assistance is not an entanglement measure according to 
some definitions, it is a very informative quantity regard- 
ing entanglement. 

Having a system of two qubits, one can also use con- 
currence instead of entanglement in Eq. (JÜJ) , yielding the 
definition of concurrence of assistance: 

C*as S ist(É>) = max^í>iC(|*i) 

i 

so that J2Pi\^i){^i\ =Q- ( 10 ) 

i 

The advantage of this quantity is, that it can be easily 
calculated for two qubits. As it is shown in [ïl|, it is 
simply 



C assist {g) = tr y/ ^fgg^fg = ^ A», 



(11) 



c.f. Eq. JB}. Notc that this quantity is essentially a fi- 
delity between the physical density matrix g and the ma- 
trix g, which is physical for separable states only. 

Thanks to the formula in Eq. (|l·lfl , concurrence of as- 
sistance is not only an informative quantity, but it is as 
feasible as concurrence itself in the case of qubit pairs. 



III. GRAPH STATES REVISITED 

In this Section we briefly review the properties of the 
Ising dynamics for spin-1/2 partides without magnètic 
field, which are known from Refs. [|J QJ. We will talk 
about spins in this context, and the eigenstates will 
represent the computational basis: |0) = |f), |1) = 
Consider a set of spins, with pairwise interactions be- 
tween them: 



H = - 



dx) 



(k,l) 



(12) 
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where the summation (k, l) goes over those spins which 
interact with each othcr. (Hence the name graph states 
for the states to be considered here: the geometry can 
be envisaged as a graph, where the vertices are the spins, 
and the edges represent pairwise Ising interactions.) As 
the summands in Eq. (|12|l commute, the time evolution 
can be written as a product of two-spin unitàries 



where 



Ú(t) = e~ i6r = JJ U k ,i{r), 
(k,i) 



(13) 



(14) 



Here r stands for the scaled time measured in arbitrary 
units. 

First we study the time instant r = ^: one may di- 
rectly verify that 



U ktl 



i + tó-W 



d*) 



(15) 



The evolution operators without a time argument will 
denote those for r = \ in what follows. These describe 
conditional phase gates in a suitably chosen basis. Let 
us assume that the system is initially in a state |e m ) of 
the computational basis, a common eigenvector of all the 



(7^ \p 

Un Cr 



~-n,m \^m/ i 



±1. 



(16) 



The state Ü \ e m ) will be an eigenvector of the following 
complete set of commuting observables: 



(17) 



with the same eigenvalues as the e n , m -s in Eq. I|16(l . The 
operators K n in Eqs. (|17|) depend on the geometry of the 
graph. They can be evaluated simply by utilizing the 
following relations: 



(x) 



Ukicr 



kJ 



(m^k,l). 



(18) 



which can be verified directly by substituting Eq. I|15|) 
into Eq. Q17JI. The joint eigenstates of these operators 
are termed as graph states [7j. It can be shown that 
many of the so arising states corresponding to different 
graphs are local unitary equivalent. 

As an example, consider a ring of N spins with pairwise 
Ising interaction. In this case 



K. 



(ring) _ _±.(x) 
l — "l-l 



d*) 



dx) 
'l+V 



(19) 



where the arithmetics in the indices is understood in the 
modulo N sense. The common eigenstates of these com- 
muting variables are termed as cluster states, and they 



were introduced in Ref. ,6j, although in a different basis. 
They are suitable as an entangled resource for one-way 
quantum computers 

Note that U(tt) = — 1 in general. Specially for a ring 
topology, U(n/2) = — i holds too. This means that the 
evolution is periòdic: at such time instants the initial 
state appears again, which is a computational basis state. 
Thus the Ising dynamics without magnètic field produces 
oscillations between the computational basis state and 
a graph (or in some of the cases, cluster) state. The 
achieved graph state is selected by the initial basis state. 

To obtain a more complete picture on the whole pro- 
cess of the entanglement oscillations, we plot the tempo- 
ral behavior of the entanglement quantities in Fig. ^ for 
the ring topology. In the figure we observe that the con- 
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FIG. 1: (Color online.) Overlap with the initial state and 
entanglement measures for the first two qubits, during the 
entanglement oscillations for five spins in a ring, generated 
by the Ising Hamiltonian without magnètic field in 11I2L In 
the initial state all spins are up, thus in state |0) if we consider 
qubits. The plotted quantities are dimensionless. 

currence of assistance of the qubit pair is almost equal 
to the square root of one-tangle of one of the constituent 
spins. We will show later in this paper that the square 
root of one tangle is an upper bound for concurrence of 
assistance. Thus for the states in argument, the entan- 
glement of a subsystem with the rest of the system can be 
indeed "focused" to a pair of qubits via suitably chosen 
measurement on the rest of the system. This is obvious 
for the cluster states, but it appears to hold for the most 
of the time evolution. 

The dynamical entanglement behavior of the systems 
in argument can be controlled by the appropriate choicc 
of the initial state. Consider for instance the following 
polarizcd initial state: 



I*a(í = 0)) 



fc=l 



|o>* 



II) 



(20) 

The "A" index reflects that all the spins are rotated from 
the z direction in the same way. This state can be pre- 
pared by a simultaneous one-qubit rotation, which is 
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available even in optical lattice systems. If 9 = Itt (l 
being integer), we obtain the graph state periodically, 
while for 8 = Itt/ 2 the state is stationary, thus no en- 
tanglement will be generated. Between these vàlues, the 
cntanglement measured by one-tangle or concurrence of 
assistance is a monotonous and continuous function of 6 
for all vàlues of time. Thus by varying this parameter 
of the initial state, one can control the amount of the 
generated entanglement. 

From the above discussion we find that Ising dynamics 
without magnètic field has the following properties from 
the point of view of entanglement generation: 

1. The generated bipartite entanglement is always 
small. 

2. In the case of the cluster states one can project 
the state with certainty to a maximally entangled 
pair of two spins by a measurement on the othcrs. 
Moreover, requircd measurement is a local one. 

3. All the states of the computational basis are peri- 
odically transferred into states which have proper- 
ties 1-2. 

4. One can control the amount of the dynamically gen- 
erated entanglement by a parameter of the initial 
state, which can be altered by the same local rota- 
tion applied on all the spins. 

During our investigations we will check which of these 
properties may arise under different couplings, initial 
states and topologies. 



IV. TWO PROPERTIES OF CONCURRENCE 
OF ASSISTANCE 

In this Section we present two properties of concur- 
rence of assistance for multi-qubit systems. 

Our first proposition formulates an upper bound of 
concurrence of assistance. 

Proposition 1 For an arbitrary state of two qubits A 
and B, square root of the one-tangle of either qubits 
serves as an upper bound for concurrence of assistance, 
i.e.: 



thus 



QA = tr B qab = y^Pfcgfc, 

k 



(24) 



due to the linearity of the partial trace. Substituting 
Eq. I|24(l into the definition in Eq. we obtain 



T A = 2 



\ 



det (^^PkQk^j , 



(25) 



while according to the definition in Eq. (|10fl . 

CÏÉ st = 2 Vdet( M ), (26) 



where we have exploited the fact that for pure states 



C7(|&» = 2 v / dit 



Qk ■ 



(27) 



Substituting Eqs. I|25(l and (|26J) into the statement of 
the Proposition in inequality l|21|l . what we have to show 
is that 



^Tdetfe 



kQk, 



de« YJ 



PkQk 



(28) 



This is a consequence of the recursive application of the 
inequality I|A1|I . which is proven in Appendix|H QED. 

Intuitively, in the spirit of the considerations concern- 
ing lower bound of localizable entanglement in Ref. ^| ; 
we can claim that a local measurement on the ancillary 
systems of a purification of qab cannot create additional 
cntanglement between the spin A and the rest of the sys- 
tem A, as such a measurement is an operation on the 
complementary system. Thus, by choosing the optimal 
measurement we can, at best, concentrate all of the orig- 
inally available entanglement (VTa) into the entangle- 
ment between the qubits A and B. 

The appearance of the one-tangle in the context of con- 
currence of assistance suggests that there might be some 
relation with the CKW inequalities, and this is the case 
indeed. Nevertheless, it is simple to prové the following: 



ÍA > C AB 



(21) 



Proof: Consider the ensemble realization of the state qab 
of the qubits A,B 



QAB = ^Pk |&) (& 



(22) 



which provides the maximum in Eq. (|10|l . and use the 
notation 



Proposition 2 For a system of three qubits A,B,C in a 
pure state, 

Cab = CIT* and C AC = C%&« (29) 

implies that the Coffmann-Kundu- Wootters inequalities 
in Eq. (jSJ are saturated, thus 



C\b + Cac — Ta 



(30) 



Qk = tr B |£ fc ) (£ fc | 



(23) 



holds 
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This immediately follows from the same derivation as 
in Ref. by exploiting the fact that the matrices Ri 
of Eq. © for subsystems AB and AC have rank one 
due to the conditions of the proposition. (C.f. Eqs. (JüJ) 
and JÏÏJO- 

Proposition [3 relates the direct and measurement as- 
sisted approach to bipartite entanglement in multipartite 
systems. The question remains open, of course, whether 
it is true for more parties, too. 

As already pointed out in Section Ulli for the graph 
states themselves y/TJ — C^ st = 1, and besides \/Ta ~ 
C*| s J> lst holds throughout the whole time evolution gen- 
erated by Ising couplings. According to Proposition 
it is correct to call such states as those with maximal 
concurrence of assistance. Meanwhile Cab <C C^ 51 , 
which suggests that CKW inequalities are far from be- 
ing saturated, which is indeed the case. The generated 
entanglement is essentially multipartite, but it can be 
converted to bipartite via a measurement. On the other 
hand, if CKW inequalities are saturated, then we can 
expect concurrence of assistance being below the square- 
root of one-tangle. Besides, the question naturally arises, 
whether it is possible to dynamically create entanglement 
oscillations in spin systems which saturate CKW inequal- 
ities instead. 



V. CONTROLLED GENERATION OF 
CONCURRENCE AND CONCURRENCE OF 
ASSISTANCE 

Now we turn our attention to spin-1/2 systems as those 
naturally realize multi-qubit systems. We assign the 



eigenstates as the computational basis states as |0) = ||), 
|1) = 14)- We will use the qubit notation for simplicity. 

We have seen in Section IIIII that certain states with 
maximal concurrence of assistance can be generated in 
dynamical oscillations, and the control over the avail- 
able entanglement is realized by the altering of the initial 
state. This control requires a simultaneous operation on 
all the spins, and as for bipartite entanglement, it effects 
the entanglement available via assistive measurements 
only, as concurrence itself takes low vàlues throughout 
the evolution. First we consider whether it is possible to 
control the concurrence itself too, and if it is possible to 
control the evolution by varying a singlo spin only. 

Consider first a system of N + 1 spins with XY cou- 
plings: 



fr — \ " a( x ) ü.( x ) i ~(y) /oi\ 
Hxy — - 2^ a i a j +<J i a j ' i ól J 

<i,j> 



in a star topology: spin is the middle one, while spins 1 
to N are the outer ones, each coupled to the central one. 
Even though the summands of the Hamiltonian do not 
commute, the eigenvalues and eigenvectors can be calcu- 
lated. One would expect that the state of the middle spin 
can control the entanglement behavior, as the interaction 
of the outer spins is mediated by this one. Indeed, if one 
considers the initial state where only the middle spin is 
rotated, the others point upwards, i.e. they are in the 
state |0): 



|*m(* = 0)) = I cos I ) |0) + sin I ~ ) |1) ) <g> I \Q) k , (32) 



fe=i 



the time evolution, as shown in Appendix[Bl reads 
|¥ M (t)) = cos (§) f|0) ® ® |0) fc 

\ k— 1 



+ sin (|) f cos(2VÏVí) |1) <g> i |0) fc - i sin(2VÏVí) |0) <g> V |0, . . . 0, l u . . .) J . 

V k=1 ^ N i=i ) 



(33) 



The rotation of the central spin indeed controls the en- 
tanglement behavior of the system: for 6 — no entan- 
glement is created, while for 9 = ir the maximal entangle- 
ment oscillation will appear. The state is a superposition 
of a product and an entangled state depending on 9, thus 
this parameter controls the available entanglement con- 
tinuously. 

These entanglement oscillations are different than 
those in case of Ising couplings. As shown in AppendixIUl 



concurrence is equal to concurrence of assistance in the 
case of any superposition of the computational basis 
states with all spins up and one down. This means 
that in the states arising throughout this evolution mea- 
surements do not facilitate "focusing" entanglement onto 
two spins. Besides, it has been proven in Ref. 3 that 
these states saturate CKW inequalities in Eq. JSJ, thus 
the bipartite entanglement present in the states is maxi- 
mal. This scheme provides a dynamical way of preparing 
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multipartite states with maximal bipartite entanglement, 
which is controlled by thc initial state of one spin. In ad- 
dition, it illustrates that Proposition [2 works for more 
than two subsystems, which is shown exactly in this spe- 
cific case. Note that at certain times the central spin gets 
disentangled from the outer ring, which is meanwhile in 
a state with highest pairwise concurrence possible. Such 
a maximally entangled state is reached for the whole sys- 
tem, too, at different times, see also in Fig. |2^a). 

In Fig. [21 we present the behavior of concurrence and 
square root of one tangles for a ring topology, and for an 
outer spin in a state different from the others, as an illus- 
tration. Here we consider the initial state producing the 
maximal entanglement, that is, one spin is considered to 
point downwards, while all the others point upwards. An 
analytical solution similar to that in Appcndix iBl would 
be feasiblc too, but more cncrgy cigcnstates have nonzero 
weights in the initial state. Of course the functions are 
not equal for all the spins in such case, but their behavior 
is similar to the star topology. According to Appendixlül 
concurrence is equal to concurrence of assistance, and of 
course CKW inequalitics are saturated. 

From the above discussion one might conclude that the 
XY couplings "prefer" to generate pure bipartite entan- 
glement. This is however not the case. In order to ex- 
amine this issue, we have plotted the behavior of entan- 
glement quantities for an XY-coupled star configuration 
with the initial state in Eq. (|2Ufl. that is, the polarized 
state arising as a product of all the spins in the same 
state which is a superposition of |0) and |1). It appears 
that in this case concurrence between two outer spins is 
heavily suppressed, but concurrence of assistance takes 
rather high vàlues for certain initial states. Moreover, 
concurrence of assistance is very close to the square-root 
of one-tangle, just as in the case of the Ising couplings. 
Thus XY couplings can, if the initial state is suitably 
chosen, produce states with a high amount of bipartite 
entanglement available via assistive measurements. No- 
tice however, that the square-root of one-tangle is higher 
than concurrence of assistance, thus there is also some 
multipartite entanglement present in the system which 
cannot be accessed by assistive measurements. 

Consider now Ising interactions, and ask whcthcr it 
is sufficient to rotate just one spin in order to control 
the amount of available entanglement, e.g. disable en- 
tanglement oscillations. For thc rotation of an outer 
spin in the star configuration or the ring topology we 
have found that entanglement cannot be completely sup- 
pressed. However, if we rotate the central spin in a star 
topology, it is possible to control entanglement behavior. 
This is illustrated in Fig. 21 Similarly to the case of initial 
state of <PÜ|) , concurrence of assistance is almost equal to 
the square root of one-tangle, while concurrence itself is 
close to zero. 

It is important to note that the possible high value of 
concurrence of assistance appears to have nothing to do 
with the bipartite nature of the couplings. In order to see 
this, consider a ring of spins with the "weird" threepartite 



couplings 

£we ir d = -£4-\4^Sr (34) 

k 

The temporal behavior of concurrence of assistance and 
square-root of one-tangle for neighbors is shown in Fig.0 
Concurrence of assistance apparently reaches its upper 
limit showing that threepartite interaction can also gen- 
erate maximal focusable bipartite entanglement. 

In this Section we have shown that it is possible to gen- 
erate entanglement oscillations not only between product 
and graph (or cluster) states, but also between product 
states, and states with maximal possible bipartite entan- 
glement, and control this entanglement behavior by the 
initial state. 



VI. ENTANGLED BASES IN THE PRESENCE 
OF A MAGNÈTIC FIELD 

In Section lïïïl we have seen that in the absence of mag- 
nètic field the Ising couplings induce such dynamics that 
all the states of the computational basis evolve into graph 
states periodically. In the Heisenberg picture we may in- 
tèrpret this so that the product of the operators 
evolves to such a joint observable, which has an eigenba- 
sis formed fully by graph states. One of the key features 
of such states is that they can be projected onto a max- 
imally entangled state of any pair of selected spins by 
a von Neumann measurement on the rest of the spins. 
We show here that this property is preserved, moreover 
enhanced if the magnètic field is present. 

First we consider the Ising Hamiltonian with a mag- 
nètic field pointing towards a direction characterized by 
the angle 4>: 

(k,l) k 

Thus we have two free parameters characterizing the 
magnètic field, its magnitude B and direction <f>. Note 
that the rotation of the magnètic field is equivalent to a 
rotation of the initial state in this case. 

In particular, we are interested in the temporal be- 
havior of the concurrence of assistance C ass i s t for certain 
pairs of spins. Therefore we calculate the time evolution 
of all the states |e<) of the computational basis: 

|eí(B ) t))=exp(-iH Wng t)|e i ) ) i = 1...2 N , (36) 

Then we can evaluate the average 

C*assïst(B, t) = Ç Cassist ^ B > *)» ' ( 3? ) 

i 

and also the Standard deviation 

*c^ ist (B,t) = ^CÏ~ t -C—? (38) 
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FIG. 2: (Color online.) Concurrence and one-tangle for spins coupled by XY interactions in the absence of magnètic field. In 
Figs. a)-d) 6+1 spins are ordered into a star topology, while in e)-f) a ring of 6 spins is considered. In the initial state all 
spins are up, except for one, which is down. In a)-b) the central spin while in c)-d) an outer spin is flipped to point upwards. 
Figures on the left display concurrences of qubit pairs, those on the right display square roots of one-tangles as a function of 
time. Legend: c: the central central spin, f: an outer spin which is flipped initially, Ofe: an outer spin which is the fc-th neighbor 
of the initially flipped one. Time is measured in arbitrary units, the other quantities are dimensionless. The figure is obtained 
from exact numerical diagonalization and direct calculations. 



of concurrence of assistance over the computational ba- 
sis states as initial states. The deviation is informative 
regarding the deviation of the quantity from the average 
for the different initial states. 

A typical result of the calculation is plotted in Fig. [Ü] 
For B = the expected entanglement oscillations are 
present. If the magnètic field is nonzero, the system does 
not tend to return to the initial product states. Mag- 
nètic field resolves many of the the high degeneracies of 
the Ising Hamiltonian, and the eigenvalues become in- 
commensurable. Therefore, even though the evolution of 
the system will be almost periòdic according to the quan- 
tum recurrence theorem [37f , the reasonable approximate 
recurrences occur after an extremely long time. 

For B ^ 0, the ensemble average of concurrence of as- 
sistance appears to be rather strictly close to one for quite 
long time intervals, while its Standard deviation is low. 
The deviation can be further suppressed by the suitable 
choice of magnètic field. This behavior of concurrence of 
assistance is very similar to that in Fig.[ü]also for different 
chosen pair of qubits, for qubit pairs of a ring topology, 
and also for different computationally feasible numbcr 
of qubits. From this we can conclude that the elements 
of the computational basis are transformed into states 
which can be projected into nearly maximally entangled 
states of chosen two spins via a von Neumann measure- 
ment on the rest of the spins. Otherwise speaking, Ising 



couplings do take the products of à^ z > matrices to such 
complete set of commuting operators, whose eigenstates 
have the above mentioned property. The temporal dura- 
tion of the presence of this property is significantly en- 
hanced by the magnètic field. 

The so arising entanglement is essentially multipartite: 
the appearance of the magnètic field does not enhancc 
concurrence of the qubit pairs as it can be verified by 
pcrforming the same calculation with concurrence. Notc 
that the characteristic behavior of the entanglement as 
reflected by the Meyer-Wallach measure for the kicked 
Ising model, also in the case of the presence of a mag- 
nètic field pointing towards an arbitrary direction was 
also reported in |24j . 

Another rclcvant qucstion might be whcthcr the re- 
quired measurements are local, i.e. how much localizable 
entanglement is present. To illustrate this issue in our 
numerical framework, we have evaluated a lower bound 
for localizable entanglement by the mere consideration of 
a measurement on the computational basis. According to 
our experience, the behavior of the so available bipartite 
entanglement resembles that of concurrence of assistance, 
but takes lower vàlues. However, quite remarkable bipar- 
tite entanglement is still available, which is in most of the 
cases still higher than the limit that CKW inequalities 
would allow for, without measurements. 

Next we investigate the properties of the AF-model 
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FIG. 3: (Color online.) Comparison of rotating all spins or the central spin in the initial state of a 6+1 spin star with XY 
couplings. Fig. a) displays the temporal behavior of concurrence if the central spin is rotated, i.e. the initial state in Eq. 11M2I is 
used, while the other three figures display the evolution of concurrence, concurrence of assistance and square-root of one-tangle 
with an initial state in Eq. 12LH . that is, all spins in the same superposition of |0) and |1). All the bipartite quantities correspond 
to two outer spins, square-root of one-tangle is that of one of these. 9 stands for the dimensionless parameter of the input state. 



from the same point of view: into Eq. (|36|1 we substitute 
the Hamiltonian 



(k,i) 



E 



i±à M ~(z) -i$-a<~' 
e a k e 



(39) 



A homogeneous magnètic field parallel to the z does not 
have any effect on the entanglement behavior of the sys- 
tem, as 



X) 



d*) 



dv) 



(k,i) 



(40) 



thus the local rotations generated by can be 

taken into account after calculating the effect of the cou- 
plings. Thcrcforc we pick B = 1, and investigate the 
dependence of concurrence and concurrence of assistance 
on the direction 4> of the field. 

The quantities evaluated are again those in Eqs. I|37|) 
and l]38p. both for concurrence and concurrence of as- 
sistance. A typical result is displayed in Fig. [7| It ap- 



pears that for 4> — we obtain oscillations in the av- 
erage concurrence, too, while concurrence of assistance 
is not significantly higher than concurrence itself. The 
appropriate choice of the direction of the magnètic field 
can suppress concurrence, significantly cnhance concur- 
rence of assistance and decrease its deviation. Thus even 
though the couplings are not Ising type, at least the fca- 
ture of the Ising couplings that it produces bases with 
high concurrence of assistance can be retained. 



VII. CONCLUSIONS 

In this paper we have related the problems of maxi- 
mizing pairwise concurrence and pairwise concurrence of 
assistance in a system of múltiple qubits. We have shown 
that the square root of one tangle of a qubit is an upper 
bound for the concurrence of assistance of a qubit pair 
containing the particular qubit. We have also shown that 
for a certain set of states for which the CKW inequality 
is known to be saturated, the concurrence is cqual to the 
concurrence of assistance. This means that the bipartite 
subsystem under consideration is not correlated with the 
rest of the system via intrinsic multipartite entanglement. 



10 





FIG. 4: (Color online.) Control of entanglement generation in a system of 6+1 Ising-coupled spins in a star configuration. The 
central spin is rotated, i.e. initial state is that in Eq. i'-Vll) . the others are in the state |0). Figures a) and c) display temporal 
behavior of concurrence as a function of parameter 9 of the initial state, for a) two outer spins and b) an outer and a central 
spin. Figure b) shows the difference between square root of one tangle and concurrence of assistance for two outer spins. Figuro 
d) shows concurrence for the central and an outer spin. This quantity is zero for the outer spins. 
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FIG. 5: (Color online.) Time evolution of concurrence of 
assistance and one-tangle for the "weird" Hamiltonian in 
Eq. (|2U, f° r 6 spins. In the initial product state all spins 
point upwards. 



We have also studied the entanglement behavior of 
spin-1/2 systems modeling qubits, from this perspective. 
We have shown that in a star configuration of an XY 
coupled spins entanglement oscillations between product 



states and states with maximal bipartite entanglement 
according to CKW inequalities can be dynamically gen- 
erated. The oscillations can be controlled by rotating the 
spin which mediates the interaction, and at some points 
it gets disentangled from the rest of the outer ring, which 
is maximally entangled in the CKW sense. This maximal 
entanglement is reached for the whole system, too. We 
have shown numerically that the star topology facilitates 
the similar control of entanglement oscillations between 
product and graph states. The rotation of all the qubits 
of the initial state on the other hand leads to different be- 
havior of concurrence of assistance, as the enhancement 
of bipartite entanglement to the measurcmcnt appears. 
We have found similar behavior for different topologies 
numerically. 

According to our numerical results magnètic field can 
lead to the temporal enhancement of concurrence of as- 
sistance in the entanglement oscillations starting from 
the states of the computational basis, in the case of spins 
coupled by Ising interactions, arranged into ring or star 
topologies. Thereby a special entangled basis can be ac- 
cessed. We have found similar behavior for the case of 
XY couplings: magnètic field applied along properly cho- 
sen direction suppresses concurrence and enhances con- 
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FIG. 6: (Color online.) Average (a,c) and Standard deviation (b,d) of concurrence of assistance for a pair of outer spins of a 
star topology, taken over all the possible computational basis states as initial states. Ising Hamiltonian with a magnètic field 
as in Eq. 135H . 4+1 spins in a ring topology. In Figs. a) and b), <f> = 0, B dependence is plotted in Figs. c) and d), B — 1, 
0-dependence is plotted. Similar figures are obtained for different choice of the spin pair, and ring topologies too. 



currence of assistance. 

According to the presented results, pairwise couplings 
between spins and qubits can be used efïectively for dif- 
ferent tasks of distributing bipartite entanglement be- 
tween múltiple parties. It is also possible to control the 
dynamical behavior of entanglement by local quantum 
operations such as rotation of control qubits. Besides, 
magnètic field can be utilized to temporally enhance cer- 
tain entanglement features, or to chose between qual- 
itatively different kinds of entanglement behavior. It 
would be also interesting to investigate whether the en- 
tangled bases available in the described means are useful 
for quantum information processing tasks. 
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APPENDIX A: AN INEQUALITY 

In this Appendix we show, that for two Hermitian, 
positive semidefinite 2x2 matrices A and B, 

VdetÀ + VdetB < y^det (A + èj (Al) 

holds. 

First we remark that for the square root of a Hermitian 
positive semidefinite matrix 

det \ÍA = V det À. (A2) 
Thus we can rewrite inequality IjAlfl as 



det \ÍÀ + det < Jdet (À + b), (A3) 
or equivalently, 

(det Vi + det V~B) 2 < det (À + b) . (A4) 

Without the loss of generality we can perform calcula- 
tions on the eigenbasis of A. Thus we can introduce the 
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FIG. 7: (Color online.) Time evolution of averages (a,c) and deviations (b,d) of concurrence (a,b) and concurrence of assistance 
(c,d) for two outer spins of a star configuration of 4+1 spins coupled by the XY Hamiltonian with magnètic field in J^J. 
Parameter <j) describes the direction of the magnètic field. Similar behavior was observed for ring topologies and different choice 
of the qubit pair too. 



notation 



A = a 



D 



ai 

1-ai 

fel fe 2 

b* 2 1-fei 



where a, fe, a\ and b\ are real. Substituting Eq. 
Eq. í|A4|> . after some calculation we obtain 



det 



(À + è) - (det y/À - det V 
= a 2 b 2 (\b 2 \ 2 + (a 1 -b 1 ) 2 )>0, 



(A5) 



H) into 



(A6) 



which is always justified. 

Note that the inequality just proven is a special prop- 
erty of 2 x 2 matrices: if we replaced A and B by positive 
numbers as "'1 x 1" matrices, the direction of inequal- 
ity (|A1(I would be reverse. 



APPENDIX B: ANALYTICAL SOLUTION FOR 
THE YY-COUPLED STAR 

Here we derive the time evolution for our specific in- 
put states in an XY coupled star configuration, based on 



Refs. [H|33|]. Consider the Hamiltonian in Eq. (JÏÏÏJ for 

a star topology of N + 1 spins. Let spin be the central 
one, thus the Hamiltonian reads 

N N 
fe=l fe=l 

Introducing the joint spin operators of the outer spins 



N 1 

Ja = 7^(1 



k=l 



a = x,y,z, 



N 



J± — ^ J x Zt íJy 



(B2) 



k=l 



and the operator for the z component of the total angular 
momcntum 

1 

the following commutation relations hold: 



L z — Jz 



(B3) 



Hxy , L 2 



H XY ,r 



0. 



(B4) 



Therefore the computational basis states with equal spins 
down span invariant subspaces of the evolution, and the 
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outer spins behave collectively as one big spin. It is 
convenient to rewrite the Hamiltonian in the Jaynes- 
Cummings type form 



Hxy = - [àí +) J- 



(B5) 



where the states \j;m) are the eigenvectors for the outer 
spins, while the states |0) and |1) are the states of the 
central spin in our qubit notation. (Note that in our 
notation, |0) = ||), thus |1) = 2 |0).) 



which has the eigenvalues and eigenvectors 

\^ m ,±) = ^=(\l) \j;m)±\Q) \j;m-l)) 



V2 

± = T2y/(j + m)(j-m+l), (B6) 



We consider the possibility of the control by the rota- 
tion of the central spin, thus our initial state reads 



|tf (i = 0)) = A |0) ® |00 . . .) + Ai |1) ® |00 . . .) = A |0) <g> 
Rewriting this in the energy basis in Eq. I|B6(I we obtain 



N 




+ Ai |1) ® 


N 




3 = Y' m = 


f) 


i = Y> m = 


f) 









(B7) 



|*(t = 0)) = A |0) 



N N\ Ai , , 



Ai 
2 



J 



N iV\ , v 



N N , 

7 = — ,m = 1 

J 2 2 



(ll>® 


TV 




-|1>® 


N 


N 




í> 


3 = ^ m = 










~2 



(B8) 



Substituting the e lwt factors, where according to Eq. I|B6|> . the u-s are 0, — 2yN, +2\/N for the three summands of 
Eq. (|B8|) respectively, after some àlgebra we obtain 



|*(í = 0)) = |0) ® |0, 0, . . .) + Ai[ cos(2V Nt) |1) ® |0, 0, . . .) - i sin(2ViVí) |0) 



J 



N N 
— , m = — — 1 
2 ' 2 



(B9) 



This shows that the complex phases of Aq and A\ are irrelevant from the point of view of the entanglement properties. 
Substituting A — cos(6/2) and A\ = sin(é'/2) into Eq. ljB9|) . and calculating \j = N/2,m — N/2 — 1) by applying 
J_ on \j = N/2,m — N/2), we obtain Eq. I|33|) . the desired result. 



APPENDIX C: RELATION OF CONCURRENCE 
AND CONCURRENCE OF ASSISTANCE FOR 
STATES WITH MAXIMUM ONE SPIN DOWN 

In this appendix we show that for states of N qubits 
of the form 



N 



i*i) = E A * 



(Cl) 



fe=0 



where 



10) = |o,o,...,o), 

\k) = |0,...0,l fc ,0...>, k^0, (C2) 

concurrence equals to concurrence of assistance for any 
pairs of the qubits. 

Consider the spins k and l. Their density matrix in 



the computational basis is of the form 

/ f?oo,oo í?oo,oi 000,10 0\ 

000,01 001,01 001,10 

000,10 001,10 0io,io 

\ oj 



(C3) 



Direct calculation of concurrence and concurrence of as- 
sistance according to Eqs. iJBJ and Eq. l(TT|) yields 



a 



kl 



2|0oi, 



Ca 
t 



2-^/001,01010. 



10- 



(C4) 



Calculating the required matrix elements from Eq. HC1|1 
we find 

0oi,oi = A* k A k , 1O) io = Ai Ai, 001,10 = A* k Ai. (C5) 



Substituting Eq. l(Ü5)l into Eq. (|Ü4|) gives C u = C í 
for arbitrary k,l. 



assist 
kl 
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